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Abstract 

The kinematic induction equation of MHD is solved numerically us- 
ing a staggered mesh in the case of the normal '111' ABC flow. Basic 
integral properties such as the growth rate of magnetic energy agree 
with previous findings, but careful 3D visualizations of the topology 
of the magnetic field reveal, however, that the conclusions about the 
modes of operation of this kinematic dynamo must be revised. The two 
known windows of dynamo action at low and high magnetic Reynolds 
number, correspond to two distinct modes, both relying on the replen- 
ishing of the magnetic field near the type stagnation points in the 
flow. The results support the case for the normal ABC flow as a fast 
dynamo. 



1 Introduction 

The class of steady velocity flows called ABC flows (after Arnold, Beltrami, 
and Childress) is an example of how complex flows may be able to amplify 
weak seed magnetic fields. ABC flows have been studied in e.g. the context of 
the cyclonic o-effect by M, but the main reason that this rather special type 
of flow is interesting to dynamo theory, is that it on the one hand constitutes 
a situation so simple that it is possible to understand the details of the 
amplification process, and that it one the other hand displays a behavior 
that in some respect resembles that of astrophysical dynamos such as the 
Sun: The amplification is related to the stretching, twisting and folding of 
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magnetic field lines, and the role played by the finite diffusivity is as essential 
to the amplification process as it is believed to be for the Sun. The main 
purpose of this paper is to illustrate the principle of operation of a dynamo 
by means of visualizing the 3-D structure of magnetic field lines and their 
dynamics. Doing so we discover yet unknown properties of this type of 
dynamo. 

2 ABC Flows 

The ABC flow is a 3-D periodic, incompressible and steady flow given by 
the sum of three parameterized Beltrami waves Q : 

u = A(0, sin kx, cos kx) + B(cos ky, 0, sin ky) + C(sin kz, cos kz, 0). (1) 

Traditionally the approach is purely kinematic, so that the back-reaction in 
the equation of motion by the Lorentz force is ignored in the treatment (but 
see H]). Thus the problem reduces to one of solving the induction equation, 

— = - V x (u x B) + ??V 2 B, (2) 

where B is the magnetic field, u is the velocity field (Eq. [I]) and r] is the 
magnetic diffusivity. 

2.1 Method 

The induction equation is solved with a simplified version of the code by 
Galsgaard, Nordlund and others (e.g. Q and ||), using both uniform mag- 
netic resistivity and the quenched diffusion that is normally used in the full 
version of this MHD code: In both cases the magnetic Reynold number Re m 
increases with the size of the grid as n 2 (where the number of grid points 
is n 3 ). In this paper only results from cases with constant and uniform 
diffusion are presented. 

The induction equation is solved on a periodic finite difference staggered 
mesh with 6th order staggering operators, 5th order centering operators, 
and a 3rd order Hyman time stepping routine. 

Even though the flow is steady, in certain regions of the flow trace par- 
ticles follow chaotic paths so that any two trace particles separate expo- 
nentially with time — a property that has been shown to be essential for 
dynamo action (||). By varying the basic parameters of the flow A, B, C, 
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and the wavenumber k it turns out that a lot of both quantitatively as well 
as qualitatively different topologies are possible: A special case with ^4 = 
is the G. O. Roberts dynamo used in the context of the Earth's dynamo 
(|0]). While the normal ABC flow with A = B = C = 1 (hence also called 
the "111" flow) has 8 stagnation points, a flow with A = 5 and B = C = 2 
has no stagnation points. If all the parameters are non-zero the flow in 
general contains a mixture of chaotic and regular islands (fl). 

The normal ABC flow has two regimes of dynamo action: One at com- 
paratively low magnetic Reynolds numbers Re m = 8.5 - 17.5 (|J) and one 
starting at Re m ps 27 (g(| and Q). 

To initialize the mode in the high Re m regime one may use a weak 
uniform magnetic seed field pointing in any direction. By flux conservation 
the uniform part remains constant, and thus rapidly becomes insignificant 
relative to the growing mode. To initialize the growing mode in the low 
Re m window, one may use the following periodic initial condition (see ]|] 
and p|), 

B = (s'm(kz) — cos(%), sin(kx) — cos(kz), sin(fey) — cos(A;x)), (3) 

which is an eigenmode of Eq. [2] in the case of zero diffusivity r\ = and is 
related to the rate of strain matrix dui/dxj of the ABC flow. 

2.2 Topology 

The normal "111" ABC flow has 8 stagnation points and 8 points where the 
velocity is maximum (see Figure |l|). A stagnation point is a point where the 
fluid velocity vanishes, and since the vorticity is aligned with the velocity 
vectors for Beltrami waves ((l3|]), the vorticity also vanishes at stagnation 
points. There are two different types of stagnation points for the normal 
ABC flow: 

• a type stagnation points where stream lines are diverging along an 
axis through the stagnation point and converging in the plane perpen- 
dicular to the axis. 

• (3 type stagnation point where stream lines are converging along the 
axis and diverging in the plane. 

The stream lines of the flow have a three-fold symmetry in the converging 
(diverging) planes through the a ((3) type stagnation points (see Figure |]). 
The three- fold symmetric 'leaves' of the converging/diverging stream lines 
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are separated by convergence/divergence connect to other stagnation points: 
While the separator lines and leaves of diverging stream lines of a f3 point 
connect to three a type stagnation points the reverse is true for an a point. 

In itself it is not so important that there happens to be 8 stagnation 
points in the normal ABC flow. After all, many interesting flows do not 
possess stagnation points (such e.g. the flow, see [14] and turbulent flows in 
general), and in any case a stagnation point may be removed by a simple 
translation of the coordinate system. Rather it is the stretching ability of 
the flow that is relevant for the dynamo action. For a divergence free flow 
V • u = such as the ABC flow the induction equation may be written in 
the form 

<9B 

— + u • VB = B • Vu + ??V 2 B (4) 

Apart from diffusive effects the stretching term BVu implies a linear growth 
of the logarithm of the magnetic field magnitude. 

In the high degree of symmetry of the normal ABC flow the stagna- 
tion points coincide with local extrema of the stretching rate. Hence the 
two types of stagnation points may be seen as convenient markers of these 
regions. 



3 Results 

Below the results of experiments with normal ABC flows are divided into re- 
sults concerning the growth of e.g. magnetic energy and exponential growth 
rates and results related to the formation and evolution of magnetic struc- 
tures. 



3.1 Magnetic energy 

As shown in Figure the initially weak seed field is indeed amplified in an 
exponential manner, and the growth increases with Re m . For intermediate 
Re m (in the second window), an oscillating behavior is associated with the 
growth (as also seen by e.g. p| , ][| and fL5|), see Figure The oscillation 
may be understood as a direct consequence of the spatial periodicity of the 
ABC flow (fU). The period of the oscillation increases with increasing Re m 
until a transition to a non-oscillating regime occurs at a Re m of about 200 
(|16|). Figure || shows examples from all three Re m regimes. 



As mentioned by [12] and |8[ it is possible to rig the initial condition of 



the seed magnetic field such that only one mode is present in the calculations; 
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the special initial condition for which that is possible is given by Eq. [3|. 

In that case the magnetic field is not amplified in the second window as 
mentioned by [ [12|| . |8| used this initial condition together with the so called 
'flux conjecture' to try to deduce the limiting growth rate of the ABC flow. 
They found exponential growth of the flux in a selected region. The growth 
rate derived does not, however, agree with the asymptotic growth rates 
found by [jl2| and fll6| . One should indeed not expect to be able to recover 
the growth rate of a completely different (exponential growing) mode by 
studying the stretching of field lines in a another (secularly decaying) mode. 

Figure || shows that the mode in the case of the initial condition given 
by Eq. || is an exponentially decaying mode: The decay is associated with 
an oscillation with a period about 10 times smaller than that of the growing 
mode. Because of numerical round-off errors the amplitude of the growing 
mode is not identically equal to zero in the initial condition, and eventu- 
ally its inevitable growth and the decay of the initialized mode results in a 
transition from decay to growth of the total magnetic energy (see Figure ^) . 



[12 1 found no growing solution and they concluded that "something odd is 



going on 



3.2 Flux cigars 

When the induction equation is evolved from a weak, uniform seed field, flux 
"cigars" rapidly arise at four of the eight stagnation points — the a-points 
(see Figure ||) . These are the regions in the flow where the magnetic field is 
most rapidly advected by a converging flow in two of the three dimensions. 
The flux cigars are aligned along the axis of divergence through the a type 
stagnation points and point directly to the (5 type stagnation points. The 
regions around the j3 type stagnation points are not unstable to flux cigar 
generation since the flow there is diverging in two directions. However, the 
converging flow along the axis from the a points creates weaker, sheet like 
flux concentrations on both sides of the plane of divergence (see Figure |8|) . 

The four flux cigars seen in Figure || correspond to the fastest growing 
eigenmode in the diffusion-less case r\ = 0. Initially several modes may be 
present but eventually the one with the largest exponential growth rate will 
be the only one that remains in the solution. 

In our experiments we do indeed find that a second set of flux cigars 
is formed next to the 4 primary cigars, but we conclude that rather than 
obscuring the physics, they are absolutely essential. These secondary cigars 
have the opposite polarity of the neighboring cigars, and they are connected 
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by reconnecting field lines (see Figure |5|) . 

The secondary flux cigars arise after the initial transient phase. These 
cigars begin to form at the a point at the separator line next to the pri- 
mary cigar. As the secondary cigar forms the primary cigar moves slightly 
away from the center of the stagnation point. In the simulation within the 
oscillatory regime, as the energy increases at the beginning of a 'cycle', the 
secondary cigar increases in size and field strength until the two cigars be- 
come equal in size midway through the cycle (Figure |6|) . The primary cigar 
now becomes smaller and actually vanishes at the end of the cycle, at which 
time the growth of the energy slows down, and even reverses briefly. At that 
point the flux cigar that was previously the secondary cigar moves into the 
center of the stagnation point (see the evolution in Figure |6|) . 

3.3 A dynamical amplification process 

The dynamical picture of an amplification cycle may be sketched as follows 
in the intermediate- Re m oscillating regime (see e.g. the case of Re m = 40 in 
Figure g): 

In the beginning of a cycle, at one of the local minima in the magnetic 
energy, there is only one flux cigar at each a point (see Figure |6|) . Gradually 
field lines of the opposite polarity pile up at the separator line close to the 
a point where a secondary cigar will eventually form. These field lines come 
from the plane of divergence of a neighboring (3 point and are transported 
along two sets of converging stream lines on opposite sides of the separator 
line. The stream lines come from two leaves of diverging stream lines at the 
[3 point separated by the separator line: Where the stream lines reach the a 
point they have twisted so that the field lines they carry are parallel to the 
axis of divergence through the a point. At the other side of the a point one 
of the leaves of diverging stream lines from a (3 point supply field lines of 
the opposite polarity to the primary flux cigar. While these two flux cigars 
with opposite polarity sit at the stagnation point and receive field lines, the 
field lines that they are made of may reconnect because of their opposite 
polarities (see in Figure ||), and the non- vanishing magnetic diffusivity. This 
reconnection process is essential to the operation of the dynamo. If there 
were no reconnections of field lines, the magnetic field near the (3 points 
could not be replenished. It is replenished by tight field line 'hooks' that 
are released to move out along the axis of divergence through the a point 
towards the plane of the (3 point. The field lines of the legs of these hooks 
that come from the primary cigar are stretched out into triangular shapes 
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by the three-fold symmetric diverging stream lines in the plane of a (3 point 
(see Figure 0). The field lines in this leg of the hook go through the plane 
of divergence of the (3 point and connect the two primary cigars on each 
side of the plane. The part of the field line that makes up the bases of the 
triangle line up perpendicular to the separator lines and are added to the 
two flux cigars at the a, points at the end of these lines. The corners of the 
triangle are also added to a points but through the diverging leaves of the 
(3 point, that become the converging leaves of three other a points. In the 
plane of the (3 points an intricate folding of field lines takes place: Figure || 
shows isosurfaces of weak magnetic field and field lines near the plane of a 
(3 point. The plane actually constitutes a region of discontinuity where the 
field lines change direction from above to below. 

Towards the end of a cycle, the primary cigar begins to decay, because 
the supply of flux to the primary cigar is less than that to the secondary 
cigar. At the end of one cycle the secondary cigar is the only remaining one 
and it moves to the center of the stagnation point. 

The magnetic energy grows through most of the cycle but the main 
growth takes place while the two cigars have about equal sizes. This is the 
time at which the most rapid reconnection take place and thus the time 
where the largest amount of magnetic flux is released down along the axis 
of divergence through the a points. 

The cycle discussed above is actually only a half cycle: A full cycle 
requires that the flux cigars return to the original polarity and this is only 
achieved after two of the above cycles. In this respect the ABC dynamo is 
similar to the solar dynamo where the full cycle of the dynamo takes two 
11 year migrations of the sunspot pattern before the polarity returns to the 
original. 

For large Re m (in the non-oscillatory regime) the amplification process 
is similar to what was described above for the lower Re m case. However, 
for large Re m , there are no oscillation associated with the amplification 
process: The double cigars always consist of one large strong cigar and one 
smaller and weaker one with the opposite polarity. The secondary cigars 
never becomes stronger than the primary but the mechanism that drives 
the amplification process remains the same. 

In the case of the single-cigar mode in the low Re m window, the mode 
pulsates with a high frequency, corresponding roughly to the time of trans- 
port of field in a logical circle. Magnetic field of the opposite polarity to 
that of a pre-existing cigar is transported towards the a point by the three 
'leaves' of converging flows. The three leaves do not provide equal amounts 
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of flux but because the diffusion is relatively large this does not result in the 
formation of a discrete cigar of opposite polarity. Rather, the existing cigar 
is surrounded from all sides by magnetic flux of opposite polarity (cf. Figure 
P), its field is cancelled, and a flux cigar of opposite polarity is formed. As 
a result of the polarity inversion, the polarity of the magnetic flux ejected 
from the (3 point region is reversed, and the cycle repeats. 

An analysis of work and dissipation shows that the average work done 
against the Lorentz force is positive throughout the cycle but that the aver- 
age Joule dissipation is smaller than the average work in the growing phase 
of the cycle and larger than the average work in the decaying phase of the 
cycle. Averaged over a cycle, the work is slightly larger than the dissipation, 
and hence the mode energy grows. 

4 Discussion and conclusions 

The dynamo action in the two windows of Re m correspond to two distinct 
modes. In both cases the replenishing of the field near the (3 points is 
crucial for the operation of the dynamo. In these regions a complicated 
folding of field lines takes place, that is accompanied by a discontinuity of 
the directions of field lines across the plane. The latter result would have 
been hard to discover without the use of modern visualization technics that 
allow direct "eye-contact" when data browsing. 

Certain properties of the magnetic transport are nearly invariant as Re m 

is increased: The size of the regions where diffusion is important become 

_i 

smaller as Re m 2 but reconnection still takes place and the field in the crucial 
(3 regions continues to be replenished. In the bulk of the flow where the 
stretching takes place the decrease of the diffusivity is unimportant since 
the field lines there are not influenced by diffusion: They tend to obtain 
a certain alignment with the flow topology given by the stretching that is 
an invariant property and hence the linear rate of increase of In B remains 
nearly the same. It thus appears very unlikely that the double-cigar mode 
should go away in the limit of infinite Re m and that the normal 111 ABC 
flow would not be a fast dynamo. 
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Figure 1: A view showing the structure of the normal ABC flow: Stagna- 
tion points are marked with purple isosurfaces, regions of maximum speed 
with blue isosurfaces, and stream lines are white. Near the center of the 
view is a (3 type stagnation point in a plane of three-fold symmetry with 
diverging stream lines. In the upper right corner is an a type stagnation 
point connecting to the (5 type point via a heteroclinic orbit; along this line 
the velocity increases from zero to a maximum half-way, decreasing again 
towards zero at the f3 point. 

Figure 2: Four panel showing the magnetic energy as function of time for 
the cases: Re m = 12, 40, 200, and 1600. 

Figure 3: The total magnetic energy Em as a function of time on a logarith- 
mic scale for two experiments (moderate Re m ) with uniform initial condition 
(thick line) and the initial condition given by Eq. ^ (thin line). 

Figure 4: A snapshot of magnetic field strength isosurfaces at a high value 
(blue). Also shown are stagnation points (green) and regions of maximum 
speed (purple). A blue flux cigar is at a green (a type) stagnation point in 
the upper corner (just above the center of the picture). This cigar points 
towards a purple region of maximum speed, a green stagnation point {(3 
type), yet another purple maximum speed region, and in the bottom corner 
is seen the opposite (periodic) end of the flux cigar. 

Figure 5: A snapshot from a high Re m experiments showing the structure of 
the double flux cigars (transparent purple) located at an a type stagnation 
point, also shown are magnetic field lines in white. 

Figure 6: Four snapshots (only part of the box is shown) of an experiment 
with a moderate Re m showing the evolution of the double flux cigars (yel- 
low) . Also shown are converging stream lines in white and the diverging axis 
through the a point (in red). The four pictures correspond to four instants 
during the oscillation in the magnetic energy, at fractions of 0, 0.14, 0.87 
and 1 of the cycle. 
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Figure 7: A snapshot from a high Re m experiment showing double flux 
cigars (purple) and field lines. The black and red field lines are connected in 
a 'hook' between the double cigars in the closer corner. The field lines are 
stretched out in a triangle by the diverging (blue) stream lines in the plane 
of the (3 point at the center of the view. 



Figure 8: A snapshot of weak magnetic field isosurfaces and field lines around 
a (3 type stagnation point. The vertical field lines near the top of the picture 
come from the double cigars at an a point. 



Figure 9: A snapshot from an experiments with Re m =12, and the initial 
condition for the magnetic field given by Eq. [|| showing magnetic field lines 
(two sets with opposite polarity) and isosurfaces, at a time near a local 
minimum of magnetic energy. The flux cigar in the center is surrounded by 
magnetic field of opposite polarity, and Joule dissipation is high between the 
two polarities. 
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